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We report the first observation of X{3872){K~^7t~) and evi¬ 
dence for X{3872){K'^tt^). The product of branching fractions for 

the former decay mode is measured to be B{B^ X(3872)(A+7r“)) x 
B{X{3872) —)• J/'ip'K^TT~) = (7.9 ± 1.3(stat.) ± 0.4(syst.)) x 10“® and 
also hnd that B^ —)• X(3872)iO*(892)° does not dominate the B^ 
A(3872)77’''7r“ decay mode in contrast to other charmonium states like 
rjj'. The product of branching fractions for the latter decay mode is mea¬ 
sured to be ^ X(3872)(A07r+)) x B{X{3872) Jfi/jTi+Ti-) = 

(10.6 ± 3.0(stat.) ± 0.9(syst.)) x 10“®. This study is based on the full and 
final data sample of 711 fb“^ (772 x 10^BB pairs) collected at the T(4S') 
resonance with the Belle detector at the KEKB collider. 
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1 Introduction 


Belle Collaboration discovered the X(3872) state [T] in the exclusive reconstruction 
of —>■ X(3872)(—)■ J/'^7r+7r“)i7+ [2] about more than a decade ago. Currently, we 

know precisely its mass (3871.69±0.17) MeV/c^ [3], have a stringent limit on its width 
(less than 1.2 MeV at 90% confidence level) |1] along with dehnitive assignment 
of [5]. It has been observed to decay to the following final states: J/'07 [S], 
|I], [8], J/'07r+7r“ [T] and D^D*^ (SKTO]. Till now, X(3872) has been 

observed and studied in two body B meson decays. This is the hrst time, we have 
observed X(3872) in three body B decay and estimated the product of branching 
fractions using full and hnal Belle data set to understand its mysterious nature. In 
this analysis, we also did the comparison of this exotic state “X(3872)” with ordinary 
charmonium states by considering tjj' as calibration sample. 

We present study of X(3872) production via the B^ —)■ X(3872)i7+7r“ and B^ 
X(3872)i7^7r+ decay modes, where the X(3872) decays to The study is 

based on 711 fb“^ of data containing 772 x 10® BB events collected with the Belle 
detector m at the KEKB e’''e asymmetric-energy collider [12] operating at the 
T(4S') resonance. 

2 Selection criterion 

To hnd the reconstruction efficiencies, signal Monte Carlo (MC) samples are gen¬ 
erated for each decay mode using EvtGen [13] and radiative effects are taken into 
account using the PHOTOS [H] package. The detector response is simulated using 
Geant3 ng. The selection criteria is same for signal MC events, background MC 
events and data events for calibration sample (having 'ip') and for concerned decay 
modes (having X(3872)) except the difference of range as both {'ip' and 

X(3872)) are further reconstructed from J/'ipTT^'K~. 

We reconstruct J/ip mesons in the decay channel (£ = e or p) and include 
bremsstrahlung photons that are within 50 mrad of either the or e~ tracks [here¬ 
inafter denoted as e+e“( 7 )]. The invariant mass of the J/'ip candidate is required 
to satisfy 3.00 GeV/c^ < < 3.13 GeV/c^ or 3.06 GeV/c^ < M^+^- < 

3.13 GeV/c^ (with a distinct lower value accounting for the residual bremsstrahlung 
in the electron mode). The J/ip candidate is then combined with a 7r’''7r“ pair to 
form an X(3872) ('^') candidate whose mass must satisfy 3.82 GeV/c^ < < 

3.92 GeV/c^ (3.64 GeV/c^ < < 3.74 GeV/c^). The dipion mass must also 

satisfy — {mj/^ + 0.2 GeV/c^), where mj/^ is nominal mass. This 

criterion corresponds to > 575 (389) MeV/c^ for the X(3872) (-0') mass re¬ 

gion and it reduces signihcantly the combinatorial background [1] with an advantage 
of flattening the background distribution in To suppress the background 
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from continuum events, we require R 2 < 0.4, where R 2 is the ratio of the second- to 
zeroth-order Fox-Wolfram moments [16] . 

To reconstruct neutral (charged) B meson candidate, a K^'k~ {Kg7i~^) candidate 
is further combined with the X(3872) for concerned decay mode and with -0' for the 
study of calibration sample. B candidates are selected using two kinematic variables: 
the energy difference AE = — i^beam and the beam-energy constrained mass 

Mbc = (F'beam ~ Pb ) /c^, where i^beam is the beam energy and E'^ and are 
the energy and magnitude of momentum, respectively, of the candidate i?-meson, 
all calculated in the e’''e“ center-of-mass (CM) frame. More details regarding the 
selection criteria can be found in Ref. na. 


3 Signal Extracton 

To extract the signal yield of R —> X(3872)(—)■ J/'iljn~^7r~)K7T, we perform a two- 
dimensional (2D) unbinned extended maximum likelihood £t to the AE and 
distributions. The 2D probability distribution function (PDF) is a product of the 
individual one-dimensional PDFs, as no sizable correlation is found. 

In order to study backgrounds, we use a large Monte Carlo sample of R —)■ J/tf^X 
events, corresponds to 100 times the integrated luminosity of the data sample. Based 
upon above study we hnd that few backgrounds are peaking in the distribution 

(nonpeaking in the AE distribution) and vice versa. The remaining backgrounds are 
combinatorial in nature and are flat in both distributions. 

For the signal, the AE dimension parametrization is done by the sum of a Crystal 
Ball [IB] and a Gaussian function while the distribution is modeled using 

the sum of two Gaussians having a common mean. The mean and resolution of 
AE and are hxed for the X(3872) mass region from signal MG samples 

after being rescaled from the results of the R° —?• decay mode. Further, 

we correct the mean of a Gaussian function for the distribution because of 

difference between the decay dynamics of ^l)' and X(3872). The tail parameters are 
hxed according to the signal MG simulation. The peaking components can be divided 
into two categories: the one peaking in but non-peaking in AE that comes 

from the R —>■ X(3872)X' decays where the X(3872) decays in J/'ijjn~^7r~ [here X' 
can be any particle], and the other peaking in AE but non-peaking in which 

comes from a R with the same hnal state where Jj ti~ is not from a X(3872). 
The peaking background in AE is found to have the same resolution as 

that of the signal, so the PDF is chosen to be the same as the signal PDF, while the 
non-peaking background in the other dimension is parameterized with a hrst-order 
Ghebyshev polynomial. For the combinatorial background in both dimensions, a hrst- 
order Ghebyshev polynomial is used. The hts are hrst validated on full simulated 
experiments and toy MG studies and no signihcant bias is seen. Fig. [T] (top) shows 
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the signal-enhanced projection plots for the —)■ X{3872){K~^7r~) decay mode. The 
result of the fit and branching fractions derived are listed in Table [T] We observe 
a clear signal for X{3872)K^7;~ with 116 ± 19 signal events corresponding 

to a significance (including systematic uncertainties related to the signal yield as 
mentioned in Table [I]) of 7.0 standard deviations (a), and measure the product of 
branching fractions to be B{B^ —> X{3872)K^7i~) x B{X{3872) —>• = 

(7.9 ± 1.3(stat.) ± 0.4(syst.)) x 10“®. The above £t is validated on the calibration 




Figure 1: Projections of (a) AE distribution with 3.860 GeV/c^ < < 

3.881 GeV/c^, (b) distribution with —11 MeV < AE < 8 MeV, (c) AE 

distribution with 3.675 GeV/c^ < < 3.695 GeV/c^, and (d) distribu¬ 

tion with —11 MeV < AE < 8 MeV for B^ X(3872)(^ J/'ipTT^Ti~)K^Ti~ decay 
mode (top) and the B^ t J/t/’7r+7r“)iF+7r“ decay mode (bottom). The curves 

show the signal [red long-dashed] and the background components [ black dashed-dot 
for the component peaking in but non-peaking in AE^ green dashed for the 

one peaking in AE but non-peaking in and magenta long dashed-dot for 

combinatorial background] as well as the overall fit [blue solid]. 


mode B^ —)■ 'ip'K^TT~. Fig. [T] (bottom) shows the signal-enhanced projection plots 
for the B^ decay mode. We measure the branching fraction to be 

B{B^ 'ip'K^TT~) = (5.79±0.14(stat.)) x lO”"^, consistent with an independent Belle 
result based on a Dalitz-plot analysis |19] . 

Further, to determine the contribution of the 77* (892) and other intermediate 
states, we perform a 2D £t to AE and in each bin of M^-k (100-MeV wide 
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Table 1: Signal yield (Y) from the fit, weighted efficiency (e) after PID correction, 
significance (S) and measnred B for X{3872)K~^tt~ and X(3872)iP°7r+. 

The first (second) nncertainty represents a statistical (systematic) contribntion. 


Decay Mode 

Yield (Y) 

^(%) 

E (<t) 

B[B X{3872)K7i)x 
B{X i3872) J/Xtt+tt-) 

B^ X{3872)K+'k- 

116 ± 19 

15.99 

7.0 

(7.9 ± 1.3 ±0.4) X lO-*^ 

5+ ^ X(3872)A:07r+ 

35 ±10 

10.31 

3.7 

(10.6 ±3.0 ±0.9) X 10-6 


bins of Mktt in the range [0.62,1.42] GeV/c^) for X (3872) mass region, which provides 
a background-subtracted Mk-k signal distribution. All parameters of the signal PDFs 
for and XE distributions are fixed from the previous 2D fit to all events. 

Then we perform a binned minimum fit to the Mk-k distribution using K*{892)^ 
and (iF^7r“)NR components, which are histogram PDFs obtained from MC samples. 
Note that the B^ X(3872)iF2*(1430)° decay is kinematically suppressed. The re¬ 
sulting fit result is shown in Fig. [2][a). We obtain 38 ± 14 (82 ± 21) signal events 
for the —)■ Y(3872)iF*(892)° (5° —)■ X(3872)(iF+7r“)NR) decay mode. This corre¬ 
sponds to a 3.0(7 significance (including systematic uncertainties related to the signal 
yield) for the 5° ^ Y(3872)(—)■ J/V’7r+7r“)iF*(892)° decay mode, and a product 
of branching fractions of B{B^ —> Y(3872)iF*(892)°) x ;B(X(3872) ^ J/'ipTT^'K~) = 
(4.0 ± 1.5(stat.) ± 0.3(syst.)) x 10“®. The ratio of branching fractions is: 

B{B^ X(3872)A:*(892)°) x B{K*{892)^ K+tt-) 

B{B^ X{3872)K+7r-) (1) 

= 0.34 ± 0.09(stat.) ± 0.02(syst.). 

In the above ratio, all systematic uncertainties cancel except those from the PDF 
model, fit bias and efficiency variation over the Dalitz distribution. 

The same procedure is also applied to the B^ —> 'ijj'K^7T~ mode. With the suf¬ 
ficient yield, we use 51-MeV wide bins of Mk-k in the range [0.600, 1.569] GeV/c^. 
We perform a binned minimum fit to the obtained Mk-w signal distribution again 
to extract the contributions of the Ktt non-resonant and resonant components. For 
this purpose, we use histogram PDFs obtained from MG samples of several possible 
components of the (iP+vr”) system: iP*(892)‘^, 7^2(1430)° and non-resonant K^'k~ 
((iP+7r“)NR). The fit result is shown in Fig. [2][b). The iP*(892)° component domi¬ 
nates and we measure B{B^ —)■ ijj' K*{892Y) = (5.88 ± 0.18(stat.)) x 10“^, which is 
consistent with the world average [5]. 

In contrast to B^ —)■ (where the ratio of branching fractions is 0.68 ± 

0.01(stat.)), B^ —)■ X(3872)7P*(892)° is not dominating in the > X(3872)7P+7r“ 
decay mode. 
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Figure 2: Fit to the background-subtracted distribution: (a) for the —)■ 

X(3872)(iF+7r“) decay mode, the curves show the —)■ X(3872)iF*(892)° [red long- 
dashed], B^ —)■ X(3872)(iF+7r“)NR [green dot-dashed], as well as the overall £t [blue 
solid], (b) for the B^ decay mode, the curves show the B^ —)■ '^'iF*(892)° 

[red long-dashed], B^ —>• [green dot-dashed], B^ (1430)° [ma¬ 

genta dashed] as well as the overall fit [blue solid]. 


We also investigate the decays B^ X(3872)(—)■ J/'07r+7r“)(X°7r+). We per¬ 
form a 2D £t to Ai7 and as before. The projections of the 2D £t for B^ 

X(3872)(—)■ J/'ijj7r^7r~){K^7r~^) in the signal-enhanced regions are shown in Figs. El^a) 
and (b). We find 35 ± 10 events for the B~^ —> X(3872)(—)■ J/'07r+7r“)(X°7r+) 
decay mode, corresponding to a 3.7 ct signihcance (including systematic uncertain¬ 
ties). The product of branching fractions is B{B'^ —)■ X(3872)X°7r+) x ;B(X(3872) 
J/'ip'K^TT~) = (10.6 ± 3.0(stat.) ±0.9(syst.)) x 10“°. The above £t is validated for the 
rj)' mass region. The projections of the 2D £t for B^ —)■ J /in 

the signal-enhanced regions are shown in Figs. [31(c) and (d). The branching fraction 
for B^ —>• '0'(^ J/'07r’''7r“)(X°7r+) is (6.00 ± 0.28(stat.)) x 10“^, while the world 
average of this quantity is (5.88 ± 0.34) x lO”"^. 

Systematic uncertainties are summarized in Table |2l All systematic uncertainties 
are added in quadrature to give total systematic uncertainty of 5.4%, 8.0%, 7.0% for 
fiO ^ X(3872)X+7r-, B+ X(3872)X07r+ and 5° ^ X(3872)X*(892)°, respec¬ 
tively. 

In summary, we report the hrst observation of the X(3872) in the decay 
X(3872)X+7r“, X(3872) ^ J. The result for the X(3872), where 5° ^ 
X(3872)X*(892)° does not dominate the 5° —)■ X(3872)(X+7r“) decay, is in marked 
contrast to the case. We have checked for a structure in the X(3872)7r and 
X(3872)X invariant masses and found no evident peaks. We measure B{B^ 
X(3872)(X+7r-)) x S(X(3872) ^ J/V'tt+tt-) = (7.9 ± 1.3(stat.) ± 0.4(syst.)) x 10“° 
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Figure 3: Projections of (a) AE distribution with 3.859 GeV/c^ < < 

3.882 GeV/c^ and (b) distribution with —11 MeV < AE < 9 MeV, (c) AE 

distribution with 3.675 GeV/c^ < < 3.695 GeV/c^, (d) distribution 

with —11 MeV < AE < 9 MeV for the X(3872)(—)■ J/'iP'k^tt~)K^'k^ decay 

mode (top) and for the J/'^7r+7r“)iF^7r=*= decay mode (bottom). Golor 

representation is same as that of neutral mode. 


Source 

X(3872) 

X(3872) 

Source 

X(3872)K:*(892)0 


K+tt- 

K^7r+ 

Lepton ID 

Kaon ID 

Pion ID 

PDF modeling 
Tracking efficiency 

Secondary B 

MG statistics 

Fit bias 

3.4 

Lepton ID 

Kaon ID 

Pion ID 

PDF modeling 
Tracking efficiency 
Kq reconstruction 
Nbb 

Secondary B 
Efficiency 

Fit bias 

3.4 

1.1 

2.5 
+1.8 
-1.3 

2.1 

1^4 

0.4 

0.6 

0.6 

3.4 

3:2 

+4.2 

-4.9 

2.5 

0.7 

1.4 

0.4 

1.0 

3.1 

1.1 

2.6 

+1.5 

-1.4 

2.1 

1.4 

0.4 

0.2 

4.6 

Total 

5.4 

8.0 

Total 

7.0 


(a) (b) 


Table 2: Summary of the systematic uncertainties in percent (a) used for 2D fit. (b) 
used for the M^-k background-subtracted fit in B^ —)■ X(3872)iF+7r“. 


and B{B+ X{m2)K^TT+) x B(X(3872) J/^/^vr+vr-) = (10.6 ± 3.0(stat.) ± 

0.9(syst.)) X 10“®. 
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